Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  982 1889 (N62W47 L= 33)   430     83        0 ( 948/1855-1018/1915)   3.4 W37-W55 N60-N62
      2  597 1685 (N61E13 L=333)   902    188        2 ( 507/1647- 683/1719)  -2.4 E26-E01 N58-N62
      3 1664 1584 (N17W63 L= 50)  3198    302        9 (1537/1530-1783/1633) -31.8 W54-W77 N10-N23
      4 1162 1544 (N31W25 L= 12)   426    122        0 (1101/1516-1217/1584)   0.0 W20-W30 N30-N32
      5 1605 1520 (N16W53 L= 40)   345     45        1 (1569/1509-1637/1536)  -9.5 W50-W56 N14-N17
      6 1758 1434 (N06W62 L= 48)  2311    274        6 (1696/1340-1820/1517)  65.6 W57-W69 N00-N11
      7  689 1450 (N40E10 L=337)   567    127        4 ( 622/1425- 754/1472) -21.8 E14-E04 N36-N41
      8 1251 1445 (N23W26 L= 12)   304     93        0 (1202/1429-1304/1464) -15.9 W22-W29 N21-N24
      9  292 1350 (N44E49 L=297)   978    146        2 ( 253/1307- 345/1415) -42.0 E56-E42 N39-N48
     10 1241 1301 (N16W20 L=  6)   232     80        0 (1229/1258-1254/1345)  53.1 W18-W21 N13-N17
     11 1835 1283 (S04W62 L= 49)  1461    119        3 (1778/1244-1882/1324)   0.0 W54-W70 S05-N00
     12  824 1224 (N22E06 L=340)   560    229        0 ( 770/1149- 896/1255)  20.6 E10-E02 N18-N24
     13 1677 1155 (S05W43 L= 29)   645    146        0 (1653/1082-1715/1220)  36.9 W39-W47 S08-S02
     14  875 1107 (N15E06 L=340)   324     99        1 ( 850/1055- 899/1155)  38.7 E08-E03 N12-N16
     15 1130  615 (S19E05 L=342)   665    211        0 (1072/ 517-1170/ 703)  36.0 E11-W00 S24-S14
     16  583  579 (S08E38 L=308)   896    189        0 ( 557/ 497- 610/ 686)  45.0 E43-E33 S12-S02
     17 1702  538 (S45W37 L= 23)  1462    187        2 (1616/ 484-1783/ 594)   0.0 W26-W48 S46-S41
     18  765  246 (S34E46 L=301)  2171    297        3 ( 639/ 166- 868/ 326)  10.7 E69-E32 S38-S30
  

    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