Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1174 1821 (N52W23 L= 64)   494    102        1 (1152/1786-1206/1850)  36.9 W19-W27 N48-N54
      2 1627 1432 (N19W44 L= 85)   648    185        5 (1575/1382-1672/1461)  -7.1 W40-W48 N16-N21
      3  799 1359 (N20E12 L= 28)   637    158        1 ( 732/1310- 888/1396) -31.0 E16-E06 N16-N22
      4 1302 1301 (N13W19 L= 59)   201     62        0 (1268/1285-1344/1313)   0.0 W16-W21 N12-N13
      5  358 1076 (N06E43 L=357)   498    140        2 ( 336/1010- 381/1153) -76.0 E44-E41 N02-N10
      6  265 1027 (N04E52 L=349)  1360    243        3 ( 214/ 906- 300/1130) -69.0 E56-E49 S03-N10
      7  128  975 (N02E69 L=331)  4365    554       11 (  74/ 873- 179/1080) -90.0 E80-E61 S03-N08
      8 1752 1019 (S06W49 L= 89)   269     43        0 (1737/ 990-1770/1049)  53.1 W47-W50 S08-S04
      9  796  362 (S43E25 L= 15)   725    153        1 ( 728/ 339- 867/ 401)   9.5 E30-E18 S44-S40
  

    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