Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1577 1525 (N33W46 L= 40)   603    118        0 (1528/1495-1630/1557) -35.5 W41-W49 N30-N35
      2  735 1333 (N24E19 L=336)  4031   1123        8 ( 328/1072-1030/1550) -27.3 E58-W00 N07-N39
      3 1023 1464 (N32W02 L=356)   248     57        0 ( 994/1447-1060/1478)  14.0 W00-W04 N31-N33
      4 1442 1412 (N27W31 L= 26)   613    204        2 (1357/1359-1530/1438) -24.4 W25-W36 N22-N28
      5 1667 1225 (N13W46 L= 40)   964    234        1 (1582/1139-1756/1295) -37.6 W39-W52 N07-N18
      6  318 1171 (N14E50 L=304)   863    177        3 ( 298/1090- 341/1270) -79.7 E51-E47 N09-N20
      7  150  928 (S02E70 L=284)   732    186        0 ( 119/ 909- 190/ 945)   8.1 E77-E63 S02-N00
      8 1385  857 (S08W22 L= 16)   476    132        2 (1330/ 817-1440/ 899) -32.0 W18-W26 S10-S04
      9  879  685 (S17E11 L=344)  1110    379        2 ( 739/ 557- 979/ 831)  43.4 E21-E03 S24-S07
     10 1537  759 (S15W33 L= 28)   278     98        2 (1503/ 744-1561/ 780)  26.6 W30-W35 S15-S13
     11 1821  641 (S26W62 L= 57)   770    204        1 (1780/ 624-1860/ 664)   0.0 W57-W67 S26-S23
     12 1061  495 (S30W01 L=355)  1823    440        6 ( 941/ 365-1243/ 601)  37.7 E09-W13 S40-S23
  

    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