Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1034 1819 (N57W35 L= 98)   399     52        0 (1003/1799-1061/1841)  10.3 W29-W40 N56-N58
      2 1437 1776 (N39W63 L=127)   488     57        0 (1415/1754-1454/1803)  10.3 W58-W69 N37-N40
      3 1464 1512 (N24W42 L=106)   412     95        0 (1421/1484-1507/1543) -51.3 W40-W44 N21-N26
      4 1093 1472 (N31W16 L= 80)   436     55        2 (1065/1439-1125/1505)  63.4 W14-W18 N29-N33
      5 1474 1203 (N07W31 L= 95)  1089    456        1 (1445/1012-1521/1402)  62.4 W25-W38 S04-N17
      6 1652 1304 (N07W47 L=111)  1248    182        3 (1583/1245-1727/1365) -66.0 W42-W52 N02-N11
      7  384 1315 (N37E41 L= 23)   271     44        0 ( 364/1296- 406/1333) -31.0 E43-E38 N35-N38
      8 1552 1301 (N10W39 L=103)   249     48        0 (1539/1275-1568/1328)  45.0 W37-W40 N08-N11
      9  931 1101 (N14E04 L= 60)   631    279        1 ( 853/1034-1004/1158) -41.2 E08-W00 N08-N16
     10 1889 1013 (S16W62 L=126)  5752    341        8 (1826/ 863-1950/1151)  41.2 W52-W75 S24-S09
     11 1817 1019 (S13W53 L=117)   612    114        1 (1782/ 973-1836/1071)  59.0 W49-W55 S16-S10
     12  441 1012 (N18E37 L= 26)  1664    401        3 ( 293/ 952- 608/1058) -24.6 E50-E26 N11-N22
     13  895 1022 (N10E07 L= 56)   230     85        0 ( 853/ 995- 932/1048) -51.3 E09-E05 N07-N12
     14 1779  902 (S19W47 L=111)  1030    157        2 (1737/ 852-1817/ 964) -58.0 W43-W52 S22-S14
     15 1280  722 (S16W08 L= 72)   837    257        0 (1194/ 635-1345/ 820)  31.6 W01-W14 S19-S11
     16 1577  517 (S39W27 L= 91)  6803    867        9 (1183/ 299-1834/ 709)   8.8 E08-W58 S44-S32
     17 1089  619 (S18E05 L= 58)   372    116        0 (1048/ 597-1153/ 653) -33.7 E07-E01 S19-S14
     18  532  534 (S12E43 L= 21)   566    110        1 ( 491/ 495- 577/ 568) -56.3 E45-E40 S15-S09
     19 1374  209 (S62W05 L= 68)   527     64        0 (1352/ 187-1402/ 237) -59.0 W02-W08 S67-S57
  

    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