Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1419 1750 (N38W53 L= 67)   497     50        2 (1397/1715-1441/1779)  26.6 W48-W56 N35-N39
      2  633 1708 (N56E17 L=357)  1780    144        4 ( 574/1648- 694/1755) -19.7 E24-E10 N50-N60
      3 1578 1717 (N31W65 L= 80)  1036    108        1 (1528/1696-1641/1733) -12.1 W58-W73 N29-N33
      4  945 1606 (N39W08 L= 23)  1399    307        0 ( 831/1536-1085/1650) -40.9 W00-W15 N31-N44
      5  363 1615 (N55E53 L=321)   558     80        1 ( 338/1584- 399/1642)   5.2 E59-E47 N51-N56
      6 1069 1406 (N23W11 L= 25)   474     80        1 (1030/1376-1111/1426)   0.0 W08-W13 N21-N25
      7  363 1312 (N32E41 L=333)  1438    314        3 ( 267/1207- 461/1392) -33.0 E53-E32 N23-N38
      8 1211 1305 (N14W17 L= 31)   595    127        0 (1139/1272-1270/1343) -45.0 W13-W19 N11-N17
      9  205 1101 (N21E58 L=317)  2843    351        7 (  94/ 997- 336/1224) -27.3 E83-E44 N13-N30
     10  428 1097 (N17E37 L=338)   403     90        2 ( 405/1066- 470/1131) -90.0 E38-E33 N15-N19
     11  532 1047 (N12E30 L=345)   317    122        2 ( 508/1013- 571/1092) -68.2 E31-E27 N09-N14
     12 1731  913 (S19W45 L= 59)   388    121        0 (1669/ 889-1783/ 951) -26.6 W39-W49 S20-S15
     13  480  814 (S00E37 L=337)   451     70        0 ( 440/ 790- 513/ 834) -31.0 E40-E35 S02-N01
     14 1330  635 (S27W11 L= 26)   426    114        1 (1288/ 587-1361/ 682)  51.3 W07-W13 S29-S24
     15  495  597 (S13E43 L=331)   281    130        0 ( 460/ 543- 520/ 656)  29.1 E48-E39 S15-S10
     16  537  534 (S18E43 L=332)   263     35        1 ( 516/ 515- 564/ 551) -56.3 E44-E40 S19-S16
     17  338  464 (S19E66 L=308)  2395    115        6 ( 293/ 408- 385/ 503) -35.5 E74-E60 S22-S15
  

    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