Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1083 1412 (N29W06 L=299)   878    179        1 (1051/1330-1117/1489) -90.0 W03-W07 N23-N33
      2 1678 1411 (N25W53 L=346)   653     96        1 (1646/1362-1706/1459) -90.0 W51-W54 N21-N27
      3  270 1414 (N30E65 L=228)   768    157        2 ( 235/1347- 322/1453)  49.4 E69-E59 N25-N32
      4  282 1347 (N25E59 L=234)   587     79        1 ( 258/1306- 306/1382)  90.0 E61-E57 N22-N27
      5 1175 1319 (N22W11 L=304)   366     76        0 (1130/1301-1218/1336)  -9.5 W08-W14 N20-N23
      6  352 1081 (N09E47 L=246)   322     95        2 ( 335/1035- 371/1121) -78.7 E48-E44 N06-N11
      7 1746  855 (S10W51 L=344)   558     60        1 (1720/ 823-1770/ 879)  36.9 W48-W53 S12-S08
      8 1029  815 (S09E00 L=293)   370     93        2 (1007/ 774-1056/ 853) -68.2 E01-W01 S11-S06
      9  841  791 (S10E12 L=281)   266     69        0 ( 825/ 756- 855/ 828)  76.0 E13-E11 S11-S07
     10 1305  748 (S14W17 L=310)   336    124        1 (1285/ 685-1320/ 801)  82.9 W15-W18 S18-S10
     11 1678  607 (S26W49 L=342)   904    142        1 (1648/ 575-1722/ 649) -35.5 W45-W54 S28-S23
     12 1669  487 (S36W54 L=347)   524     61        1 (1640/ 461-1692/ 504) -14.0 W49-W57 S37-S34
  

    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