Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  718 1806 (N44E38 L= 26)   623     44        0 ( 691/1780- 747/1830)  56.3 E40-E35 N41-N46
      2 1455 1293 (N16W24 L= 88)   299     41        0 (1424/1267-1483/1316) -26.6 W21-W25 N14-N16
      3  235 1094 (S08E56 L=  9)   746    180        2 ( 181/1044- 282/1151) -45.0 E61-E50 S10-S04
      4  453 1032 (S10E36 L= 28)   350     91        1 ( 435/ 989- 472/1086)  78.7 E37-E35 S12-S07
      5  217 1021 (S12E57 L=  7)   488     85        1 ( 201/ 984- 241/1065)  78.7 E58-E54 S14-S09
      6 1644  846 (S08W42 L=107)  1199    244        1 (1556/ 811-1727/ 891)   4.1 W35-W49 S09-S04
      7  380  786 (S26E42 L= 22)   577    173        0 ( 352/ 726- 432/ 858) -74.1 E44-E38 S29-S21
      8 1074  740 (S22W06 L= 71)   210     69        0 (1048/ 701-1095/ 777)  68.2 W05-W07 S24-S19
      9 1470  585 (S26W36 L=101)   360     65        0 (1435/ 558-1495/ 621)  68.2 W34-W37 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