Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  530 1554 (N11E47 L=  6)   355     79        1 ( 495/1529- 557/1582)  36.9 E48-E44 N08-N12
      2  672 1507 (N11E35 L= 18)   234     57        0 ( 634/1488- 701/1529)  53.1 E36-E32 N09-N13
      3 1155 1446 (N21E05 L= 48)   347    150        0 (1098/1379-1196/1506) -29.1 E09-W00 N17-N23
      4 1584 1185 (N19W28 L= 81)   735    257        1 (1559/1050-1610/1304) -60.3 W24-W32 N11-N25
      5  300 1267 (S10E54 L=359)   315     71        3 ( 279/1231- 314/1294) -53.1 E56-E52 S12-S08
      6  675 1236 (S04E26 L= 28)   203     82        0 ( 653/1191- 692/1277) -78.7 E26-E24 S06-S01
      7 1741 1017 (N15W44 L= 97)  1088    201        1 (1711/ 927-1762/1119) -70.0 W39-W47 N08-N19
      8  360  890 (S31E40 L= 13)  1315    436        3 ( 214/ 858- 508/ 934)   6.3 E54-E27 S33-S28
      9 1051  715 (S23W11 L= 64)   301     59        1 (1019/ 691-1079/ 748)  90.0 W09-W11 S24-S20
     10 1207  432 (S34W32 L= 85)  1412    180        5 (1086/ 381-1346/ 526)  60.6 W26-W36 S39-S23
     11 1502  427 (S24W51 L=105)  1228    219        1 (1390/ 390-1578/ 487)  45.0 W44-W55 S28-S18
     12 1621  381 (S21W65 L=118)   582     82        0 (1580/ 351-1662/ 417)  56.3 W62-W66 S24-S18
  

    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