Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  468 1195 (N15E37 L=204)   543    168        3 ( 442/1131- 500/1283) -60.9 E40-E34 N11-N20
      2  775 1235 (N18E16 L=226)   284    125        1 ( 702/1212- 818/1269)   7.1 E20-E12 N16-N19
      3  322 1135 (N11E49 L=192)   823    209        0 ( 300/1027- 355/1237)  85.6 E51-E46 N04-N17
      4  860  973 (N01E10 L=231)   400    115        1 ( 828/ 910- 886/1033) -66.8 E12-E08 S02-N05
      5  687  934 (S01E21 L=220)   213     96        0 ( 646/ 886- 724/ 959) -36.9 E23-E19 S03-N00
      6 1508  746 (S15W31 L=273)   304     95        0 (1488/ 697-1538/ 796)  66.8 W30-W33 S18-S11
      7  772  670 (S18E17 L=224)   849    241        3 ( 729/ 560- 814/ 780) -86.2 E19-E14 S25-S10
  

    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