Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1003 1588 (N34W03 L= 10)  4208    759        7 ( 699/1377-1288/1739) -35.3 E22-W19 N18-N49
      2 1495 1551 (N28W39 L= 46)  2265    472        0 (1341/1425-1665/1714) -55.7 W31-W49 N18-N41
      3 1566 1667 (N36W52 L= 59)   542     97        0 (1516/1646-1609/1690) -20.6 W47-W56 N34-N38
      4  426 1527 (N35E43 L=324)  1230    307        4 ( 305/1485- 564/1566) -12.3 E61-E29 N30-N39
      5  745 1501 (N30E16 L=352)   540    175        1 ( 681/1452- 815/1538) -31.0 E20-E10 N26-N33
      6  767 1356 (N20E14 L=353)   499    186        1 ( 721/1323- 824/1414) -35.5 E17-E10 N18-N24
      7  205 1326 (N24E62 L=305)  2069    165        3 ( 151/1258- 249/1399) -26.6 E80-E55 N18-N29
      8  505 1311 (N20E32 L=335)   297     83        0 ( 465/1279- 539/1335) -26.6 E35-E29 N17-N21
      9  814 1172 (N09E12 L=355)   255    162        0 ( 747/1151- 872/1219)  14.0 E15-E07 N07-N11
     10  281 1183 (N14E51 L=317)   443    112        0 ( 243/1151- 327/1211) -26.6 E54-E46 N11-N15
     11 1000  999 (S03E02 L=  6)   203     53        0 ( 970/ 973-1022/1012) -33.7 E03-W00 S04-S02
     12  915  848 (S11E08 L=359)  1159    317        3 ( 872/ 700- 988/ 981)  64.8 E12-E02 S20-S03
     13  364  730 (S14E48 L=319)   919    182        0 ( 333/ 642- 406/ 807)  42.3 E53-E42 S19-S09
     14 1116  735 (S20W04 L= 11)   508    150        0 (1073/ 664-1159/ 788)  49.4 W00-W06 S23-S16
     15  261  569 (S23E67 L=300)  2026    190        3 ( 228/ 484- 294/ 648)  54.0 E74-E60 S28-S17
     16  567  428 (S36E42 L=325)   585    103        1 ( 538/ 383- 599/ 475)  31.0 E47-E37 S38-S32
     17  433  337 (S41E67 L=300)   983     72        1 ( 390/ 318- 466/ 359)   3.4 E76-E59 S41-S39
  

    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