Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  837 1810 (N36E41 L=166)   759    193        1 ( 804/1750- 867/1859) -30.3 E47-E35 N32-N40
      2 1284 1619 (N37E03 L=205)   246     37        0 (1260/1602-1312/1635) -14.0 E04-W00 N35-N38
      3  521 1077 (S16E31 L=176)   820    202        1 ( 437/ 990- 588/1143)  61.7 E35-E27 S23-S10
      4 1825 1017 (N18W52 L=260)   568    102        0 (1800/ 969-1844/1066) -39.8 W48-W54 N15-N20
      5  996  965 (S10W00 L=208)   313    103        0 ( 954/ 921-1029/ 999) -36.9 E02-W03 S11-S08
      6  651  921 (S22E19 L=189)   283    108        0 ( 628/ 887- 694/ 975) -18.4 E21-E15 S23-S19
      7 1743  882 (N08W49 L=256)  1239    249        0 (1699/ 789-1803/ 965) -17.4 W42-W58 N04-N11
      8 1261  857 (S08W18 L=225)   594    155        2 (1238/ 786-1284/ 950) -66.0 W15-W19 S12-S03
      9 1196  794 (S14W16 L=224)   205    109        0 (1131/ 763-1230/ 819)   0.0 W11-W19 S14-S12
     10  446  689 (S42E30 L=177)   980    295        0 ( 341/ 667- 560/ 711)  19.3 E40-E20 S45-S37
     11 1033  567 (S31W15 L=223)   221    101        0 ( 987/ 539-1085/ 595)  -7.1 W11-W19 S31-S30
     12 1275  530 (S26W31 L=239)   651    244        0 (1167/ 496-1408/ 557)  32.7 W25-W39 S30-S21
  

    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