Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  647 1772 (N33E52 L=261)   527    148        2 ( 615/1734- 668/1798) -20.6 E55-E47 N30-N35
      2  379 1521 (N11E58 L=254)  5634    913       14 ( 293/1311- 496/1747) -71.6 E66-E53 S03-N27
      3  455 1499 (N11E50 L=262)   594    133        3 ( 425/1456- 483/1528) -48.2 E53-E47 N08-N13
      4  428 1093 (S13E39 L=274)   578     89        2 ( 401/1039- 463/1134)  80.5 E40-E36 S16-S10
      5  970  771 (S19W03 L=316)   503    181        1 ( 891/ 744-1072/ 799)  33.7 E01-W09 S22-S16
      6  295  687 (S42E47 L=265)  1005    290        3 ( 234/ 638- 362/ 731)  29.4 E55-E39 S46-S37
      7 1711  576 (S10W61 L= 13)  1019    200        4 (1653/ 544-1763/ 605)   8.7 W54-W67 S12-S08
      8 1392  448 (S27W42 L=354)   575    184        1 (1315/ 420-1460/ 481)  15.3 W35-W46 S28-S25
      9 1350  267 (S39W51 L=  4)  1106    276        1 (1327/ 211-1380/ 343) -31.6 W45-W58 S42-S34
  

    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