Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1750 1413 (N05W59 L=266)   395     55        2 (1721/1389-1784/1439) -31.0 W56-W62 N03-N07
      2 1440 1387 (N13W34 L=241)   317     63        0 (1404/1365-1470/1419)   0.0 W31-W37 N12-N14
      3 1655 1309 (N03W46 L=253)   398    104        0 (1626/1264-1688/1354)  20.6 W42-W50 N01-N04
      4 1181 1169 (N08W12 L=219)   743    188        3 (1125/1097-1259/1222)  -6.3 W08-W17 N04-N11
      5  154 1138 (N31E65 L=143)   799     90        0 ( 123/1105- 193/1170) -18.4 E70-E58 N28-N33
      6 1763 1122 (S10W49 L=256)   298     95        1 (1740/1079-1786/1168) -80.5 W47-W50 S13-S07
      7  401 1008 (N19E39 L=168)   397     98        1 ( 377/ 956- 424/1053)  45.0 E41-E36 N15-N20
      8  787  877 (N02E17 L=190)   408    147        1 ( 731/ 817- 830/ 929) -74.1 E18-E15 S01-N06
      9 1471  734 (S23W19 L=226)  1503    368        2 (1374/ 589-1544/ 880)  35.2 W09-W26 S29-S16
     10 1042  826 (S07E04 L=203)   452    209        0 ( 965/ 795-1107/ 872) -40.6 E07-W00 S08-S02
     11  703  705 (S05E26 L=181)   456    186        0 ( 677/ 623- 748/ 764)  63.4 E28-E22 S10-S02
     12 1878  732 (S38W60 L=267)   813     77        2 (1844/ 708-1914/ 757) -14.9 W53-W68 S39-S35
     13  978  650 (S15E12 L=195)   590    235        1 ( 924/ 567-1048/ 713)  26.6 E17-E07 S19-S12
     14  315  552 (S06E62 L=146)  4845    412        4 ( 209/ 379- 383/ 696)  78.7 E80-E54 S17-N04
     15  828  564 (S16E24 L=184)   257     12        0 ( 804/ 546- 852/ 585)  14.0 E25-E21 S17-S14
     16 1276  438 (S37E01 L=206)   368     92        0 (1249/ 399-1305/ 476)  51.3 E03-W01 S39-S34
     17 1013  427 (S29E18 L=189)   234     77        1 ( 994/ 397-1041/ 471)  80.5 E19-E16 S32-S26
     18  570  335 (S24E53 L=154)  1392    183        3 ( 499/ 296- 659/ 385)   2.9 E63-E42 S25-S22
  

    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