Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  615 1718 (N63E10 L=299)   400     74        0 ( 591/1683- 641/1748)  29.7 E13-E05 N60-N65
      2 1481 1679 (N28W55 L=  4)   343     82        0 (1444/1658-1519/1702)  -5.7 W50-W60 N27-N28
      3 1102 1665 (N40W28 L=337)   379     32        1 (1077/1641-1128/1687)   0.0 W25-W30 N38-N41
      4  962 1617 (N42W16 L=324)   244     60        1 ( 933/1599- 994/1638)   9.5 W13-W19 N40-N42
      5 1321 1498 (N24W32 L=341)   415    140        0 (1282/1439-1372/1525)  14.0 W28-W36 N20-N24
      6 1403 1440 (N18W35 L=343)   414    133        0 (1375/1381-1424/1494)  39.8 W30-W37 N15-N21
      7 1862 1173 (S10W61 L=  9)  1261    203        3 (1814/1089-1911/1264)  21.8 W53-W70 S13-S06
      8  317 1114 (N28E45 L=263)  1081    241        0 ( 232/1042- 412/1173) -35.2 E54-E37 N21-N33
      9 1637  859 (S20W33 L=341)  2073    423        7 (1524/ 694-1740/1052)  34.3 W22-W44 S27-S12
     10  198  813 (N13E64 L=245)  1056    166        3 ( 144/ 754- 254/ 869) -38.7 E70-E58 N08-N16
     11 1127  770 (S10E01 L=308)   638    137        0 (1100/ 697-1152/ 844)  63.4 E03-W02 S14-S06
     12  268  787 (N10E56 L=252)   529    111        1 ( 237/ 739- 300/ 832) -71.6 E58-E53 N06-N12
     13 1087  541 (S23E10 L=299)   303     87        0 (1064/ 512-1120/ 586) -80.5 E10-E08 S25-S19
     14 1226  460 (S32E04 L=304)   696    101        2 (1177/ 432-1285/ 484) -23.2 E07-W00 S34-S29
     15 1010  320 (S36E24 L=285)  1577    265        3 ( 893/ 276-1130/ 387)  -3.0 E33-E13 S37-S33
     16  746  159 (S41E60 L=249)   721     67        1 ( 722/ 131- 763/ 188)  17.4 E69-E53 S43-S38
  

    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