Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  537 1734 (N25E58 L=330)   822    113        1 ( 509/1686- 559/1784) -48.8 E62-E53 N21-N29
      2  421 1684 (N20E65 L=322)  1432     98        2 ( 368/1643- 464/1717)  26.6 E71-E60 N16-N22
      3 1151 1628 (N32E09 L= 19)   562    142        0 (1133/1566-1181/1700) -52.1 E12-E04 N28-N37
      4 1619 1418 (N33W29 L= 57)  1703    284        2 (1529/1297-1701/1523) -29.4 W20-W37 N26-N38
      5 1203 1453 (N22E00 L= 27)   206     70        0 (1187/1420-1228/1488) -56.3 E01-W01 N20-N24
      6 1712 1214 (N22W39 L= 67)   390     59        0 (1675/1190-1741/1237)   8.1 W35-W42 N21-N23
      7  488  949 (S23E32 L=356)  1127    226        3 ( 398/ 897- 600/ 993)  32.7 E38-E24 S28-S18
      8  896  707 (S28W01 L= 28)   330    166        0 ( 829/ 681- 982/ 728)  32.0 E03-W05 S31-S25
      9  573  671 (S39E22 L=  6)   478    104        1 ( 534/ 627- 617/ 716)  49.4 E24-E18 S42-S34
     10  395  424 (S62E51 L=337)   596     87        0 ( 361/ 403- 445/ 453)  -3.8 E56-E40 S63-S60
     11 1179  329 (S43W37 L= 65)   424    119        0 (1120/ 294-1227/ 360)   0.0 W30-W43 S44-S42
  

    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