Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  641 1788 (N32E56 L=261)  1202    162        0 ( 593/1747- 706/1824)  -3.8 E63-E46 N29-N34
      2  618 1612 (N20E43 L=273)   950    116        1 ( 589/1548- 649/1675) -66.0 E46-E40 N15-N24
      3 1415 1449 (N31W13 L=329)  2124    533        2 (1302/1215-1513/1632) -42.4 W00-W24 N19-N40
      4  294 1548 (N09E70 L=247)  3241    206        2 ( 220/1481- 356/1626)  63.4 E79-E65 N04-N14
      5  345 1475 (N05E58 L=258)   381     57        2 ( 325/1450- 366/1505) -36.9 E60-E55 N04-N07
      6  872 1398 (N13E19 L=298)   218     61        0 ( 852/1362- 884/1423) -36.9 E20-E16 N10-N14
      7 1656 1300 (N29W34 L=350)   805    136        1 (1597/1250-1706/1336) -11.3 W28-W38 N27-N30
      8 1099 1156 (N05W01 L=318)   764    231        1 (1011/1087-1149/1249)  74.7 E02-W02 N00-N11
      9 1468 1108 (N12W23 L=340)   261    122        0 (1425/1063-1487/1145) -53.1 W20-W24 N09-N13
     10 1727 1016 (N14W44 L=  1)  1614    277        3 (1643/ 925-1780/1102) -41.6 W36-W49 N09-N18
     11  629  991 (S16E22 L=295)   258     76        0 ( 605/ 957- 650/1022)  90.0 E23-E20 S18-S14
     12 1605  845 (N01W39 L=356)   743    154        2 (1530/ 816-1667/ 888)   4.8 W32-W44 N00-N03
     13 1220  746 (S15W19 L=335)   259     87        0 (1193/ 713-1244/ 788) -78.7 W17-W19 S16-S11
     14 1459  717 (S10W34 L=351)   247     19        1 (1437/ 696-1477/ 735)  90.0 W32-W35 S10-S08
     15 1289  689 (S16W24 L=341)   579    161        0 (1221/ 657-1374/ 713)  36.9 W21-W29 S19-S13
     16  252  618 (S49E59 L=258)  2420    209        4 ( 164/ 572- 335/ 686)  -3.7 E79-E45 S52-S45
     17 1658  587 (S11W54 L= 11)   560    147        0 (1617/ 550-1688/ 646)  90.0 W51-W56 S13-S07
     18 1057  362 (S42W26 L=342)   760    126        2 (1002/ 332-1118/ 397)  12.5 W21-W30 S43-S40
  

    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