Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  750 1587 (N27E32 L=228)  2433    374        7 ( 659/1436- 816/1725) -35.3 E46-E22 N18-N35
      2 1189 1706 (N47E06 L=255)   466     93        0 (1136/1693-1238/1723)  14.0 E09-E01 N45-N48
      3  567 1383 (N11E37 L=224)   368     97        1 ( 542/1346- 613/1411)  45.0 E38-E33 N08-N12
      4 1677 1290 (N28W41 L=301)   623    139        0 (1625/1241-1736/1320)  -5.2 W35-W46 N26-N29
      5 1255 1211 (N14W10 L=270)   582    170        2 (1222/1136-1286/1260) -45.0 W07-W12 N10-N17
      6 1030  856 (S11W04 L=264)   601    111        2 ( 977/ 832-1083/ 888)   0.0 W00-W07 S13-S08
      7 1083  699 (S20W10 L=271)   494    195        1 (1008/ 665-1160/ 735)  41.2 W06-W14 S22-S15
      8  444  636 (S38E35 L=225)   321     48        1 ( 419/ 612- 463/ 664) -90.0 E37-E33 S39-S35
      9  469  389 (S58E41 L=219)  1161    171        1 ( 400/ 366- 545/ 414)  10.1 E56-E27 S61-S54
     10 1004  371 (S44W16 L=276)   489    135        0 ( 940/ 352-1060/ 400)  45.0 W11-W19 S46-S39
     11 1327  353 (S36W42 L=302)   828    151        1 (1258/ 331-1386/ 382)  18.4 W35-W46 S38-S33
  

    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