Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  993 1907 (N62W50 L=199)   536     51        0 ( 976/1879-1008/1940)  10.8 W42-W63 N60-N65
      2  432 1462 (N47E33 L=116)   422     41        0 ( 411/1435- 458/1485)  26.6 E35-E30 N44-N48
      3 1029 1332 (N22W08 L=158)   187     75        0 ( 991/1316-1066/1344) -11.3 W05-W10 N20-N22
      4  851 1292 (N24E03 L=146)   233     68        1 ( 835/1251- 874/1322)  56.3 E04-E01 N22-N25
      5  337 1124 (N27E43 L=106)   772    187        1 ( 304/1044- 400/1177)  45.0 E47-E37 N22-N29
      6 1690  901 (S20W38 L=187)   524    109        3 (1654/ 845-1747/ 948)  16.7 W34-W44 S22-S18
      7 1618  679 (S32W29 L=178)   510    178        0 (1540/ 633-1685/ 729)   8.7 W22-W35 S32-S30
      8  848  568 (S17E22 L=127)   240     96        0 ( 802/ 540- 885/ 599)   8.1 E26-E19 S17-S15
      9  439  542 (S09E51 L= 98)   633     97        1 ( 421/ 485- 467/ 585)  40.6 E55-E47 S12-S06
     10 1418  495 (S38W11 L=160)   648    149        0 (1353/ 454-1478/ 541)  10.3 W05-W16 S39-S35
     11  890  304 (S35E33 L=117)   321     67        0 ( 865/ 275- 908/ 339)  23.2 E36-E29 S36-S32
     12  560  286 (S28E60 L= 89)  1579    142        3 ( 503/ 253- 626/ 336) -50.2 E64-E53 S31-S23
  

    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