Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1613 1533 (N24W56 L=246)   728    145        3 (1577/1476-1649/1596)  76.0 W52-W58 N20-N28
      2 1436 1524 (N28W41 L=230)   840    320        2 (1304/1488-1582/1579) -29.1 W32-W51 N22-N33
      3 1046 1336 (N25W08 L=197)  1133    282        2 ( 957/1232-1117/1437) -71.6 W03-W10 N17-N32
      4  893 1405 (N32E01 L=189)   232     70        0 ( 864/1371- 921/1436)  36.9 E03-W01 N30-N34
      5  617 1106 (N18E24 L=165)   321    107        0 ( 586/1064- 642/1162) -74.1 E26-E22 N15-N22
      6  833 1045 (N11E11 L=179)   325    146        1 ( 803/ 999- 897/1078)  23.2 E13-E06 N08-N13
      7  960  813 (S04E08 L=182)   451    124        0 ( 914/ 780-1017/ 851) -50.2 E09-E04 S07-S01
      8 1091  677 (S15E02 L=187)   502    144        1 (1061/ 616-1112/ 748)  58.0 E05-W00 S18-S10
      9 1282  592 (S25W08 L=198)  1173    445        0 (1146/ 471-1443/ 738)  27.6 E03-W20 S30-S18
     10 1445  578 (S30W20 L=209)   906    288        1 (1356/ 476-1532/ 665)  32.7 W12-W26 S35-S25
     11  348  534 (S14E61 L=129)   951    110        2 ( 299/ 502- 398/ 562)  -5.2 E66-E54 S15-S11
     12  823  477 (S24E24 L=165)   415     68        0 ( 792/ 444- 848/ 516)  38.7 E27-E21 S25-S21
  

    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