Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  597 1530 (N33E36 L=331)   450     92        2 ( 578/1476- 615/1581) -49.4 E38-E32 N29-N36
      2 1587 1480 (N33W43 L= 50)   832    138        2 (1519/1438-1639/1518) -21.8 W37-W47 N30-N36
      3  871 1480 (N31E13 L=354)   559    223        1 ( 787/1420- 953/1517) -17.1 E19-E06 N27-N33
      4 1318 1434 (N30W19 L= 26)   391    208        0 (1222/1388-1399/1462) -14.0 W12-W24 N26-N31
      5 1114 1199 (N14W05 L= 12)   321    104        1 (1083/1180-1177/1227)   0.0 W03-W09 N12-N15
      6 1796 1043 (N06W56 L= 62)   441    112        0 (1756/1004-1839/1067)  20.6 W51-W60 N03-N06
      7  475  933 (S05E36 L=331)   466     74        1 ( 451/ 888- 502/ 970) -78.7 E37-E33 S08-S03
      8 1689  904 (S03W46 L= 53)   786    120        2 (1638/ 867-1740/ 947) -18.4 W41-W50 S05-N00
      9 1057  800 (S11W03 L= 10)   278     98        1 (1018/ 764-1092/ 828) -38.7 W00-W05 S13-S09
     10 1815  766 (S12W61 L= 68)   461     54        2 (1795/ 743-1840/ 796) -45.0 W58-W65 S13-S09
     11 1413  770 (S12W26 L= 33)   232     63        1 (1376/ 752-1446/ 783)  11.3 W23-W28 S12-S10
     12  865  675 (S20E09 L=358)   931    223        5 ( 805/ 596- 913/ 766) -70.0 E12-E05 S24-S13
     13  621  719 (S18E26 L=341)   215    105        1 ( 589/ 696- 683/ 742)  23.2 E28-E21 S19-S16
     14  975  693 (S18E02 L=  5)   621    116        4 ( 931/ 641-1020/ 737) -51.3 E04-W01 S21-S15
     15  419  566 (S30E45 L=321)  1544    288        2 ( 359/ 504- 533/ 665) -31.6 E49-E35 S34-S23
     16 1154  486 (S32W12 L= 18)  4031    604       10 ( 954/ 390-1383/ 648)  28.2 E02-W26 S39-S19
     17  755  581 (S27E17 L=350)   790    188        1 ( 692/ 531- 834/ 630)  35.0 E21-E11 S30-S22
     18 1510  498 (S30W39 L= 46)  1587    339        3 (1381/ 432-1664/ 568)  25.5 W29-W51 S35-S24
  

    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