Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1148 1866 (N53E25 L=255)   440     38        0 (1119/1853-1177/1878)  53.1 E27-E22 N50-N54
      2  518 1727 (N23E60 L=220)   595     70        1 ( 484/1701- 538/1749)  18.4 E64-E56 N21-N24
      3  796 1594 (N20E30 L=250)   242     35        1 ( 779/1571- 815/1613) -33.7 E31-E28 N18-N21
      4 1341 1463 (N27W07 L=287)   758    151        1 (1294/1390-1390/1537) -45.0 W02-W11 N22-N30
      5 1810 1133 (N22W49 L=329)   449     39        1 (1788/1109-1836/1156) -11.3 W46-W52 N20-N23
      6 1755 1084 (N17W45 L=324)   422     90        0 (1732/1030-1772/1113) -45.0 W42-W46 N14-N19
      7  555  837 (S29E25 L=255)   216     86        0 ( 519/ 813- 594/ 867)  45.0 E27-E22 S31-S26
      8 1714  648 (S07W56 L=336)   335     36        1 (1697/ 623-1729/ 669)  90.0 W54-W58 S08-S05
      9 1513  434 (S25W53 L=333)   400    100        0 (1481/ 397-1559/ 474)  80.5 W52-W55 S27-S21
     10 1499  323 (S31W63 L=342)   432     57        1 (1469/ 296-1521/ 351)  -5.2 W57-W68 S31-S29
  

    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