Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1106 1513 (N28W11 L= 72)   391    191        0 (1015/1479-1184/1534) -26.6 W05-W15 N25-N30
      2  792 1415 (N25E11 L= 50)   822    202        1 ( 688/1366- 864/1443) -24.4 E17-E06 N21-N28
      3 1148 1321 (N16W11 L= 72)   263    105        1 (1107/1304-1208/1346) -26.6 W08-W14 N14-N17
      4  234 1254 (N22E56 L=  5)  1729    287        1 ( 156/1164- 352/1340) -25.7 E70-E43 N14-N28
      5  658 1026 (N03E22 L= 39)   242     67        0 ( 625/1009- 691/1043) -26.6 E24-E19 N02-N04
      6 1746  950 (S12W48 L=109)   611    139        3 (1714/ 890-1776/1017)  66.8 W45-W50 S15-S08
      7  380  969 (N03E42 L= 19)   317     49        0 ( 366/ 931- 397/ 999)  63.4 E43-E40 N00-N04
      8 1197  891 (S10W09 L= 70)   256    112        0 (1159/ 861-1248/ 925)  18.4 W06-W12 S12-S08
      9 1667  735 (S25W42 L=103)   971    220        1 (1640/ 640-1688/ 857)  90.0 W39-W45 S31-S17
     10 1759  667 (S30W53 L=114)   885    114        1 (1711/ 636-1797/ 711)  18.4 W47-W57 S32-S27
     11  368  491 (S26E59 L=  2)  1722    209        2 ( 324/ 385- 407/ 583)  31.8 E74-E51 S32-S19
  

    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