Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1356 1762 (N55W41 L=163)   688     96        2 (1320/1727-1393/1795)  10.6 W34-W50 N51-N57
      2 1315 1561 (N38W25 L=147)   388     54        2 (1282/1535-1346/1582) -21.8 W22-W27 N36-N39
      3  578 1424 (N30E32 L= 90)  1999    427        3 ( 426/1288- 714/1537) -27.3 E51-E20 N20-N37
      4 1068 1256 (N18W03 L=126)   437    189        0 (1028/1179-1128/1310)  53.1 W00-W07 N13-N21
      5  755  973 (N01E17 L=105)  1401    432        3 ( 664/ 817- 883/1106)  48.6 E23-E08 S08-N09
      6 1481  921 (S04W29 L=152)   249     76        0 (1466/ 887-1503/ 961) -63.4 W28-W30 S05-S01
      7  630  637 (S20E28 L= 95)  1393    338        1 ( 544/ 496- 711/ 738)  43.3 E37-E20 S30-S13
      8 1413  514 (S31W28 L=150)  3955    719        8 (1142/ 357-1611/ 726)  39.8 W07-W41 S42-S16
      9  953  671 (S18E05 L=117)   228     56        1 ( 933/ 652- 978/ 692)  33.7 E06-E03 S19-S16
     10 1019  511 (S29E01 L=121)  1908    500        6 ( 846/ 373-1201/ 670)  35.8 E15-W11 S39-S18
  

    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