Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1449 1633 (N42W14 L= 17)  1395    203        3 (1355/1575-1545/1689)   0.0 W05-W22 N39-N44
      2  722 1525 (N16E30 L=333)   523    155        0 ( 698/1452- 742/1587) -53.1 E33-E26 N11-N19
      3 1579 1522 (N38W26 L= 29)   546     82        2 (1552/1491-1605/1555)  59.0 W23-W28 N35-N40
      4 1330 1282 (N15W12 L= 15)   437    138        0 (1310/1205-1354/1351) -63.4 W09-W14 N10-N18
      5  366 1065 (S18E43 L=320)   640    159        0 ( 340/ 978- 390/1142) -66.0 E46-E40 S22-S13
      6 1870 1053 (N17W57 L= 60)   546     90        4 (1851/1018-1896/1087) -76.0 W55-W59 N14-N19
      7  882  983 (S13E07 L=356)   896    212        0 ( 788/ 887- 949/1084)  74.1 E10-E04 S19-S05
      8 1588 1044 (N08W33 L= 36)   242     65        0 (1551/1031-1618/1057)  21.8 W30-W35 N06-N08
      9 1494  823 (S07W33 L= 35)  1291    335        1 (1358/ 687-1608/ 972)  69.1 W28-W36 S17-N04
     10 1173  685 (S23W17 L= 20)   943    229        3 (1089/ 652-1276/ 748)   4.4 W10-W23 S24-S20
     11 1681  669 (S10W52 L= 55)   611    110        2 (1639/ 622-1713/ 722)  81.9 W49-W53 S13-S06
  

    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