Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1629 1536 (N22W53 L=221)   712     88        0 (1586/1498-1673/1573) -51.3 W50-W55 N19-N24
      2  604 1503 (N36E23 L=145)   884    232        0 ( 522/1439- 714/1532)   7.1 E29-E13 N33-N38
      3 1173  722 (S20W05 L=172)   467    155        0 (1114/ 669-1217/ 782)  36.9 W00-W08 S22-S16
      4  275  613 (S13E62 L=105)  1496    157        3 ( 229/ 561- 310/ 678)  81.9 E66-E57 S16-S09
      5  969  651 (S21E09 L=159)   198     76        0 ( 929/ 637-1007/ 669) -31.0 E11-E06 S22-S19
      6  333  609 (S15E56 L=112)   630    106        0 ( 314/ 560- 355/ 662)  80.5 E58-E54 S17-S11
      7  660  593 (S20E30 L=137)   217     85        0 ( 640/ 558- 682/ 635)  45.0 E32-E27 S22-S17
      8  395  385 (S30E67 L=101)   623     55        0 ( 379/ 349- 408/ 418)  39.8 E71-E62 S32-S27
  

    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