Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  899 1759 (N49E03 L= 65)  1106    312        2 ( 804/1656-1001/1847)  41.4 E10-W08 N41-N57
      2  639 1519 (N33E23 L= 45)   827    130        2 ( 590/1477- 700/1554)   0.0 E26-E18 N30-N35
      3 1407 1455 (N22W29 L= 97)   433    129        1 (1344/1422-1466/1477) -20.6 W24-W32 N19-N23
      4 1193 1356 (N17W13 L= 81)   240     76        1 (1160/1317-1236/1391) -59.0 W11-W15 N14-N19
      5  376 1313 (N22E42 L= 26)   487    130        1 ( 337/1286- 422/1350) -35.5 E45-E38 N19-N24
      6 1351 1335 (N15W23 L= 91)   287    132        1 (1287/1317-1417/1348)  -6.3 W18-W27 N13-N15
      7 1816 1285 (N08W59 L=127)   794    104        2 (1789/1233-1837/1325) -90.0 W56-W62 N05-N11
      8  231 1021 (N05E55 L= 12)  1178    273        1 ( 177/ 907- 269/1145)  47.5 E62-E50 S01-N12
      9  970  815 (S14E05 L= 63)   608    139        1 ( 933/ 757-1020/ 885)  39.8 E07-E01 S17-S09
     10  556  329 (S42E47 L= 20)   557    116        1 ( 511/ 311- 606/ 350)  -6.3 E51-E42 S43-S40
  

    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