Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1171 1799 (N51W44 L= 62)   528     73        1 (1143/1769-1213/1817)   5.2 W38-W49 N49-N52
      2  228 1264 (N36E61 L=318)   728     93        0 ( 196/1218- 261/1300) -28.6 E67-E55 N32-N38
      3  801 1038 (N13E13 L=  6)  1531    419        2 ( 730/ 923- 937/1199) -75.3 E15-E06 N03-N22
      4  338 1076 (N23E46 L=332)   510    155        0 ( 278/1045- 393/1138) -32.5 E52-E41 N20-N27
      5  401  893 (N11E43 L=335)   595    131        3 ( 358/ 855- 446/ 933) -33.7 E46-E39 N08-N13
      6 1510  844 (S15W26 L= 44)   311     89        0 (1490/ 798-1529/ 882)  45.0 W23-W27 S17-S12
      7  193  814 (N08E67 L=312)   711     88        0 ( 163/ 791- 230/ 843) -26.6 E70-E62 N06-N10
      8 1071  546 (S23E08 L= 10)   366     91        0 (1041/ 504-1101/ 588)  38.7 E11-E05 S24-S20
      9  721  324 (S31E41 L=338)   442     78        0 ( 702/ 284- 747/ 350)  56.3 E43-E37 S34-S28
     10  579  272 (S32E59 L=319)  1248    110        3 ( 537/ 217- 607/ 318)  51.3 E67-E54 S38-S28
  

    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