Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  730 1808 (N44E37 L= 26)   787     90        1 ( 700/1776- 764/1838)  53.1 E41-E33 N40-N46
      2  894 1515 (N23E14 L= 49)   278     38        0 ( 874/1486- 912/1540) -45.0 E15-E12 N21-N24
      3  388 1472 (N16E49 L= 14)   404     88        1 ( 374/1427- 407/1512) -78.7 E51-E47 N13-N18
      4 1464 1291 (N16W24 L= 88)   248     38        0 (1434/1266-1488/1312) -26.6 W22-W26 N14-N16
      5  242 1090 (S08E55 L=  8)  1574    353        3 ( 201/ 953- 287/1227)  76.0 E59-E52 S16-N00
      6  469 1032 (S10E35 L= 28)   320     89        1 ( 451/ 989- 482/1081)  68.2 E36-E34 S12-S07
      7 1653  845 (S08W43 L=106)  1313    268        1 (1570/ 806-1742/ 901)  -3.8 W35-W50 S09-S04
      8  391  782 (S26E41 L= 22)   589    180        0 ( 361/ 727- 447/ 853) -50.2 E43-E36 S29-S21
      9 1102  761 (S20W08 L= 71)   245     79        0 (1082/ 722-1124/ 806)  78.7 W06-W08 S22-S17
     10 1480  577 (S26W37 L=101)   406     60        0 (1446/ 548-1510/ 614)  68.2 W35-W38 S28-S23
  

    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