Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1291 1755 (N53W36 L= 54)   684     96        1 (1248/1727-1342/1772)   4.4 W30-W43 N50-N54
      2 1581 1603 (N37W55 L= 73)  1322    223        1 (1508/1552-1685/1656) -34.7 W48-W64 N32-N41
      3 1725 1568 (N33W72 L= 90)   741     53        2 (1699/1545-1751/1587)   3.6 W66-W82 N31-N33
      4  524 1483 (N36E38 L=341)   361     89        1 ( 498/1446- 547/1523)  71.6 E39-E36 N33-N39
      5  670 1432 (N32E24 L=354)   841    194        1 ( 613/1376- 746/1504) -36.9 E29-E17 N28-N37
      6 1522 1321 (N20W36 L= 55)   459    223        1 (1484/1265-1584/1395) -69.4 W33-W40 N15-N24
      7  772 1187 (N16E15 L=  3)  1019    467        3 ( 632/1087- 941/1259) -16.7 E24-E04 N09-N20
      8 1214 1147 (N11W13 L= 31)   240     46        1 (1181/1132-1244/1163)   0.0 W10-W14 N10-N12
      9 1193 1073 (N07W11 L= 29)   262     82        1 (1161/1043-1218/1102)  36.9 W08-W12 N05-N08
     10 1559  811 (S12W34 L= 52)   529    138        1 (1517/ 761-1618/ 852)  29.7 W31-W38 S15-S10
     11  317  510 (S27E65 L=313)   554     91        0 ( 292/ 470- 338/ 550) -90.0 E67-E62 S29-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