Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1370 1740 (N47W45 L=131)  1890    323        5 (1331/1648-1421/1819)  52.4 W39-W53 N39-N52
      2 1147 1504 (N34W16 L=102)  1064    220        0 (1067/1408-1237/1575) -65.2 W11-W19 N26-N39
      3  308 1343 (N33E51 L= 34)   360    108        0 ( 285/1320- 342/1374) -26.6 E53-E47 N31-N35
      4 1738 1250 (N08W53 L=139)  1180    308        1 (1644/1174-1834/1324) -34.5 W46-W62 N02-N13
      5 1316 1214 (N13W21 L=106)   342    115        1 (1258/1179-1359/1247)  14.0 W16-W24 N12-N14
      6  234 1080 (N16E57 L= 29)   464    171        1 ( 185/1045- 272/1106) -18.4 E62-E53 N14-N18
      7  385 1043 (N13E43 L= 43)   679    189        0 ( 324/ 977- 441/1101) -41.6 E47-E38 N09-N17
      8 1583  832 (S14W34 L=120)   393    120        1 (1543/ 808-1633/ 859) -33.7 W31-W37 S15-S11
      9  430  706 (S08E44 L= 42)   351     56        2 ( 403/ 685- 459/ 728) -11.3 E46-E41 S08-S06
     10 1585  586 (S31W34 L=120)  2344    437        5 (1403/ 494-1729/ 713)  18.4 W19-W46 S34-S24
     11  809  492 (S25E22 L= 64)   772    144        2 ( 747/ 461- 871/ 533)  16.7 E26-E16 S26-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