Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  669 1710 (N24E46 L=240)   888    113        1 ( 622/1677- 726/1736)  -5.2 E50-E39 N21-N25
      2 1338 1618 (N38W01 L=287)  1244    368        1 (1176/1545-1490/1699)   2.2 E12-W14 N36-N39
      3  367 1517 (N06E58 L=227)   444     78        1 ( 339/1478- 395/1555) -78.7 E61-E55 N03-N08
      4  766 1461 (N11E27 L=258)   733    261        1 ( 698/1378- 833/1549) -23.2 E34-E20 N07-N15
      5 1607 1247 (N23W29 L=314)   374    113        0 (1560/1220-1655/1271)   8.1 W25-W32 N21-N24
      6  135 1123 (S21E70 L=215)   781    140        1 ( 105/1057- 175/1169)  41.6 E75-E65 S25-S17
      7 1836 1138 (N25W51 L=337)   357     98        0 (1797/1109-1871/1153)  12.5 W47-W56 N22-N26
      8 1326  982 (S00W18 L=303)   244     72        0 (1293/ 963-1352/1012)  56.3 W15-W18 S02-N01
      9  869  952 (S14E07 L=279)   268     43        1 ( 845/ 925- 898/ 979) -76.0 E07-E04 S16-S12
     10 1886  845 (N11W65 L=350)  2284    176        4 (1849/ 778-1937/ 933) -48.8 W59-W78 N07-N15
     11 1228  801 (S13W18 L=303)   454    144        1 (1186/ 740-1269/ 871) -40.6 W13-W21 S15-S09
     12 1132  718 (S20W15 L=300)   338    120        0 (1065/ 700-1187/ 731)  20.6 W10-W18 S21-S18
     13 1841  705 (N02W66 L=351)   820     81        0 (1797/ 685-1881/ 727)  20.0 W61-W72 N00-N04
     14  272  655 (S48E53 L=232)   459     53        0 ( 243/ 641- 308/ 673)  10.3 E58-E47 S49-S46
     15 1038  641 (S27W12 L=298)   455    189        1 ( 949/ 605-1106/ 665)  15.3 W06-W17 S29-S25
     16 1584  515 (S16W53 L=338)  1143    221        3 (1487/ 480-1686/ 551)  34.7 W46-W60 S21-S12
     17  716  463 (S48E02 L=283)   423     56        1 ( 681/ 439- 747/ 488)   0.0 E05-W00 S49-S46
     18 1000  372 (S44W23 L=308)   434     98        1 ( 962/ 344-1049/ 405)  -6.3 W19-W28 S45-S42
     19 1502  297 (S30W65 L=351)   557     67        1 (1471/ 271-1529/ 322)   0.0 W59-W71 S31-S29
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y