Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1422 1684 (N33W47 L= 75)  1965    318        4 (1279/1617-1545/1734) -40.9 W36-W55 N26-N39
      2  507 1670 (N56E33 L=355)  2313    224        8 ( 442/1601- 572/1720)  52.1 E43-E23 N50-N60
      3  774 1554 (N39E05 L= 23)  1296    243        3 ( 685/1482- 881/1599) -47.5 E11-W01 N32-N44
      4  874 1342 (N23E03 L= 25)   470     74        1 ( 846/1312- 915/1376)  14.0 E04-W00 N21-N25
      5  340 1234 (N27E43 L=345)   486    104        3 ( 312/1189- 368/1280) -50.2 E46-E40 N24-N30
      6  995 1234 (N15W02 L= 30)   466    142        0 ( 931/1201-1053/1279) -54.5 W00-W05 N11-N18
      7  405 1134 (N20E38 L=350)   847    141        3 ( 377/1067- 446/1197)  66.8 E40-E34 N16-N23
      8  304 1156 (N23E47 L=341)   548     77        3 ( 272/1127- 338/1183)  -8.1 E50-E43 N21-N24
      9 1803 1068 (S11W53 L= 81)  1142    176        0 (1780/ 969-1831/1148)  48.4 W49-W57 S16-S06
     10  156 1015 (N17E65 L=323)  1746    236        4 (  85/ 955- 242/1112) -20.3 E81-E54 N13-N23
     11  310 1041 (N16E47 L=341)  1049    163        3 ( 258/ 992- 373/1082) -29.1 E52-E41 N12-N18
     12 1636  851 (S20W36 L= 64)   386    116        0 (1590/ 812-1695/ 881) -41.2 W32-W40 S24-S17
     13  289  761 (N00E55 L=333)   419     61        0 ( 255/ 736- 319/ 788) -31.0 E57-E51 S01-N02
     14  869  561 (S22E19 L=  9)   385    171        0 ( 820/ 490- 903/ 631)  37.9 E24-E15 S25-S18
     15 1140  564 (S28E03 L= 26)   371     95        1 (1103/ 521-1173/ 603)  38.7 E05-W00 S29-S24
     16  416  472 (S19E57 L=331)   888    113        1 ( 385/ 430- 460/ 504) -63.4 E60-E53 S22-S16
  

    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