Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1531 1589 (N21W52 L=208)  1188    189        4 (1458/1537-1624/1623) -14.0 W44-W60 N17-N23
      2 1139 1571 (N33W25 L=181)   493     76        0 (1101/1545-1180/1595)   0.0 W21-W28 N31-N34
      3 1757 1350 (N02W58 L=214)   983    154        1 (1720/1266-1788/1422)  58.0 W55-W60 S03-N05
      4  683 1395 (N36E11 L=145)   218     53        0 ( 651/1379- 714/1405) -21.8 E13-E08 N35-N37
      5  148  898 (N19E70 L= 86)   549     69        3 ( 131/ 853- 164/ 929) -59.0 E72-E66 N15-N20
      6 1304  771 (S16W09 L=165)   511    221        1 (1233/ 702-1375/ 862)  24.4 W03-W14 S18-S12
      7  842  692 (S08E19 L=137)   966    250        0 ( 766/ 591- 916/ 797) -86.2 E21-E17 S15-N00
      8 1332  559 (S30W05 L=161)   380     96        0 (1285/ 543-1382/ 576) -26.6 W02-W08 S31-S27
      9  978  369 (S31E24 L=132)   989    165        1 ( 919/ 292-1018/ 432)  76.0 E27-E20 S37-S25
  

    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