Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  249 1331 (N02E63 L= 49)   530     88        0 ( 215/1307- 275/1356)  26.6 E66-E60 N00-N03
      2  969 1292 (N13E09 L=103)   435    133        1 ( 929/1253-1020/1356) -29.7 E12-E05 N11-N15
      3  833 1117 (S00E13 L= 99)   247     56        0 ( 816/1083- 854/1148) -33.7 E14-E11 S01-N01
      4  383 1041 (S13E42 L= 70)   502    143        0 ( 364/ 959- 400/1104) -83.7 E43-E40 S18-S09
      5 1125  952 (S04W07 L=119)   349     82        1 (1097/ 904-1146/1003)  90.0 W06-W08 S06-N00
      6  945  931 (S09E03 L=109)   382    116        1 ( 917/ 872- 969/ 989) -81.9 E04-E01 S12-S05
      7  257  888 (S25E51 L= 60)   401    193        1 ( 204/ 867- 313/ 904)   4.8 E57-E45 S26-S23
      8  793  787 (S20E09 L=103)   543    221        0 ( 686/ 750- 893/ 835)  36.9 E15-E03 S24-S15
      9  794  679 (S27E07 L=105)   465    199        0 ( 714/ 623- 881/ 730)  48.0 E11-E02 S32-S22
     10 1153  680 (S19W15 L=127)   377    138        0 (1109/ 624-1210/ 723)  74.1 W13-W18 S23-S16
     11 1280  689 (S16W23 L=135)   216     60        1 (1255/ 663-1316/ 708)  45.0 W21-W24 S17-S14
     12 1348  319 (S38W44 L=156)  1021    190        1 (1270/ 301-1424/ 344)  20.0 W37-W49 S40-S34
  

    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