Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  861 1824 (N63E23 L= 13)   970    145        2 ( 795/1804- 933/1849)  -5.4 E32-E11 N60-N65
      2 1163 1728 (N53W14 L= 50)   375     40        2 (1135/1707-1192/1745) -11.3 W11-W17 N50-N54
      3  564 1471 (N32E36 L=  1)   459    138        1 ( 523/1427- 603/1539) -41.6 E41-E31 N28-N36
      4 1208 1389 (N27W13 L= 49)   801    159        1 (1141/1339-1256/1426) -29.7 W08-W15 N23-N29
      5  875 1306 (N21E10 L= 26)   369    159        0 ( 823/1276- 920/1342) -23.2 E13-E06 N19-N23
      6 1752  854 (S08W52 L= 89)   545     38        2 (1727/ 830-1777/ 881) -21.8 W49-W55 S09-S06
      7 1414  761 (S13W25 L= 62)   626    202        0 (1386/ 699-1463/ 837) -77.5 W23-W29 S17-S08
      8  847  686 (S18E12 L= 25)   789    163        2 ( 776/ 651- 924/ 718)  21.8 E16-E06 S20-S15
      9  393  522 (S31E52 L=344)  2653    342        8 ( 312/ 434- 491/ 636) -55.3 E58-E44 S37-S22
     10 1423  576 (S26W28 L= 65)  1259    328        0 (1310/ 537-1541/ 635)   6.7 W20-W37 S28-S21
  

    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