Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  569 1662 (N36E38 L= 87)  3041    537        4 ( 419/1592- 736/1742) -16.4 E57-E23 N30-N43
      2  495 1419 (N19E37 L= 88)   888    269        4 ( 438/1369- 566/1473) -15.3 E41-E30 N15-N22
      3  389 1306 (N11E43 L= 81)   784    179        2 ( 372/1221- 408/1393) -79.7 E45-E40 N06-N17
      4 1179 1230 (N09W09 L=133)   506    144        0 (1138/1163-1227/1304) -58.0 W05-W11 N04-N13
      5  711 1048 (S04E19 L=106)   212     47        0 ( 694/1022- 728/1081) -76.0 E20-E17 S05-S01
      6 1825  720 (S17W61 L=186)  1403    308        1 (1727/ 621-1899/ 844)  47.0 W53-W68 S24-S09
      7 1747  789 (S14W51 L=176)   890    163        1 (1696/ 750-1808/ 826)  29.1 W47-W56 S16-S11
      8  504  629 (S30E36 L= 89)  1408    448        2 ( 409/ 480- 645/ 784)  42.3 E46-E23 S39-S19
      9 1193  687 (S24W12 L=137)   247     60        1 (1156/ 669-1227/ 706)  26.6 W10-W14 S25-S22
     10  918  566 (S33E05 L=119)   332     45        1 ( 895/ 542- 942/ 591)  45.0 E07-E03 S34-S31
  

    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