Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  637 1765 (N27E51 L=249)   579    155        1 ( 616/1725- 659/1803) -29.1 E56-E47 N24-N29
      2  926 1749 (N34E31 L=269)   292     49        0 ( 901/1728- 952/1770) -11.3 E33-E28 N33-N35
      3  907 1367 (N09E16 L=284)  1037    299        5 ( 883/1212- 936/1527) -60.9 E21-E11 N00-N18
      4 1113 1176 (N04W01 L=301)   363    187        0 (1060/1105-1150/1216)  90.0 E02-W02 N00-N06
      5  584 1077 (S15E27 L=273)  1410    470        1 ( 351/1020- 753/1111)  15.9 E44-E16 S19-S11
      6  225 1090 (S22E56 L=245)   391     55        2 ( 200/1071- 251/1108)  11.3 E58-E53 S23-S20
      7  137  998 (S29E63 L=237)  2050    255        1 (  97/ 953- 183/1048)  40.6 E67-E59 S31-S25
      8 1679  535 (S12W59 L=360)   393    141        2 (1662/ 500-1707/ 572) -45.0 W56-W63 S14-S10
      9 1746  503 (S11W69 L= 10)   744    108        2 (1729/ 471-1770/ 533) -45.0 W66-W73 S12-S09
     10 1682  467 (S15W64 L=  5)   692    186        2 (1644/ 438-1723/ 504)   5.7 W58-W68 S16-S13
     11 1302  398 (S32W40 L=340)  1222    251        2 (1188/ 360-1424/ 462)  57.3 W33-W44 S38-S24
     12 1552  298 (S28W65 L=  5)   627    218        0 (1524/ 272-1578/ 324)   6.3 W60-W69 S28-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