Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1367 1623 (N41W33 L=331)   277     79        0 (1348/1595-1397/1643)  18.4 W30-W36 N39-N42
      2 1344 1458 (N30W26 L=323)   338    101        0 (1303/1443-1399/1471)   0.0 W22-W29 N29-N30
      3 1086 1257 (N19W05 L=303)   219     89        0 (1058/1215-1123/1295) -56.3 W03-W07 N15-N21
      4  467 1232 (N20E37 L=260)   239     73        1 ( 448/1204- 483/1260) -63.4 E38-E35 N17-N21
      5 1707  775 (S16W47 L=344)  1989    446        3 (1540/ 616-1808/ 953)  39.3 W33-W56 S25-S05
      6  847  914 (S01E12 L=286)   192     59        0 ( 815/ 899- 884/ 928)  14.0 E13-E09 S02-N00
      7  304  825 (S05E52 L=245)   648     81        2 ( 280/ 786- 334/ 868) -68.2 E54-E49 S07-S02
      8 1826  752 (S19W59 L=357)   780     96        0 (1808/ 706-1844/ 813)  74.1 W57-W61 S22-S15
      9 1870  713 (S23W66 L=  4)   860    105        1 (1847/ 664-1891/ 760)  56.3 W63-W68 S26-S20
     10  655  524 (S26E29 L=268)   287     52        1 ( 624/ 501- 682/ 546) -33.7 E31-E26 S27-S24
  

    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