Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  891 1466 (N15E19 L=238)  1259    299        2 ( 844/1339- 936/1607) -49.1 E26-E13 N07-N22
      2  314 1443 (N02E59 L=197)  1521    315        2 ( 229/1403- 410/1479)  17.4 E68-E50 N00-N05
      3 1501 1324 (N22W21 L=277)   429     86        0 (1478/1275-1519/1366) -80.5 W18-W22 N19-N25
      4 1037 1255 (N06E05 L=252)   654    138        2 ( 965/1212-1100/1308) -11.3 E10-W00 N04-N07
      5  213 1110 (S19E58 L=199)  2513    465        6 ( 159/1025- 264/1219) -70.0 E63-E51 S23-S12
      6 1496 1062 (N07W26 L=283)  1337    265        1 (1458/ 925-1541/1202) -51.8 W21-W32 N00-N14
      7  652  863 (S25E19 L=238)   859    286        1 ( 526/ 830- 784/ 914)  32.0 E27-E11 S29-S19
      8  264  785 (S38E50 L=207)  1552    435        1 ( 172/ 727- 392/ 820)  22.2 E59-E37 S42-S33
      9 1645  408 (S22W63 L=320)  1835    156        2 (1586/ 361-1694/ 473)  71.6 W59-W67 S26-S17
  

    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