Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1579 1547 (N25W48 L=312)  1001    222        3 (1530/1460-1629/1623) -90.0 W44-W51 N19-N30
      2 1121 1532 (N30W13 L=277)   734    250        1 (1029/1435-1188/1604) -61.7 W08-W15 N22-N35
      3  727 1419 (N27E15 L=249)   422    151        0 ( 673/1379- 793/1468) -45.0 E18-E10 N23-N31
      4 1916  922 (S16W69 L=333)   577     43        2 (1899/ 897-1931/ 952) -45.0 W66-W72 S17-S14
      5 1815  901 (S17W55 L=319)   389     76        0 (1796/ 870-1847/ 927)  26.6 W53-W59 S18-S14
      6 1705  623 (S33W48 L=312)  1273    252        3 (1610/ 583-1808/ 670)  -2.5 W38-W62 S35-S30
      7 1354  559 (S33W17 L=281)   554    136        0 (1324/ 508-1373/ 627)  76.0 W15-W18 S36-S28
      8 1461  515 (S38W26 L=291)   427    125        0 (1409/ 483-1506/ 560)  29.7 W22-W29 S39-S35
      9 1043  438 (S38E07 L=257)   331     75        0 (1009/ 410-1068/ 475)  31.0 E09-E04 S39-S35
     10  794  419 (S36E26 L=238)   421     72        0 ( 763/ 394- 834/ 449) -45.0 E27-E23 S37-S33
     11  864  232 (S52E31 L=233)   386     54        0 ( 833/ 219- 895/ 246)   8.1 E34-E27 S53-S50
  

    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