Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  631 1791 (N32E57 L=261)  1204    163        1 ( 581/1754- 696/1826)  -9.5 E65-E47 N28-N34
      2  601 1616 (N19E45 L=273)   930    137        0 ( 578/1546- 626/1677) -60.9 E48-E41 N14-N23
      3  296 1557 (N10E71 L=247)  2936    201        5 ( 225/1485- 354/1635)  61.4 E78-E66 N04-N15
      4 1402 1454 (N31W12 L=330)  2125    523        4 (1300/1221-1498/1632) -43.7 W00-W22 N19-N40
      5 1645 1304 (N29W33 L=350)   860    126        1 (1584/1253-1695/1344) -11.3 W27-W37 N26-N30
      6 1076 1173 (N05E01 L=317)   932    358        0 ( 951/1095-1132/1325)  78.7 E06-W01 S01-N14
      7 1452 1118 (N12W22 L=340)   266    114        0 (1414/1070-1471/1153) -45.0 W19-W23 N09-N13
      8 1712 1039 (N15W42 L=  0)  1333    222        3 (1626/ 950-1770/1109) -13.0 W35-W48 N11-N19
      9  609 1020 (S15E24 L=294)   212     81        1 ( 580/ 993- 633/1047) -14.0 E26-E22 S16-S13
     10 1588  854 (N01W38 L=356)   769    153        0 (1516/ 823-1658/ 898)   4.8 W31-W43 N00-N03
     11 1339  759 (S11W25 L=343)   218     59        1 (1310/ 737-1369/ 778)  53.1 W23-W26 S12-S08
     12 1330  706 (S14W26 L=344)   822    278        0 (1226/ 674-1466/ 739)  34.7 W20-W34 S18-S08
     13  248  619 (S49E60 L=258)  2585    212        5 ( 159/ 577- 335/ 684)  -1.6 E81-E45 S52-S45
     14 1658  602 (S10W54 L= 11)   452    116        1 (1627/ 556-1680/ 655) -63.4 W50-W55 S13-S07
     15 1292  501 (S26W33 L=351)   238     75        0 (1259/ 471-1332/ 518)   0.0 W30-W36 S27-S25
     16  442  472 (S55E34 L=284)   334     36        1 ( 425/ 451- 460/ 500)  59.0 E36-E31 S57-S52
     17 1038  367 (S43W24 L=342)   635     95        2 ( 981/ 342-1089/ 393)  14.0 W19-W28 S44-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