Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  424 1570 (N11E55 L=297)   387    104       -1 ( 404/1537- 443/1597) -23.2 E58-E51 N09-N12
      2 1540 1518 (N37W20 L= 13)   271     81        1 (1509/1502-1576/1529)  38.7 W18-W23 N34-N39
      3 1768 1363 (N35W42 L= 34)   847    140        2 (1725/1314-1808/1400) -68.2 W37-W46 N31-N38
      4  279 1167 (S15E53 L=300)   725    195        2 ( 222/1107- 326/1229)  53.1 E56-E48 S19-S10
      5  268 1091 (S20E52 L=301)  1849    665        2 ( 104/1042- 439/1164)  -4.6 E73-E36 S22-S16
      6 1569 1104 (N12W30 L= 22)   243     93        0 (1550/1061-1582/1149) -53.1 W27-W31 N09-N13
      7 1843  924 (N11W57 L= 49)  1331    305        6 (1802/ 866-1903/ 987)  56.3 W52-W63 N06-N15
      8 1454  606 (S18W38 L= 31)   757    211        2 (1437/ 518-1469/ 706) -45.0 W33-W42 S22-S13
      9 1089  554 (S32W18 L= 10)  1362    422        2 ( 881/ 490-1293/ 602)  11.7 W03-W32 S35-S27
  

    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