Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  885 1863 (N42E45 L=148)   925    124        2 ( 835/1817- 941/1893)  10.3 E52-E37 N39-N44
      2  500 1677 (N18E58 L=136)   480     76        1 ( 482/1644- 521/1709) -59.0 E60-E54 N15-N20
      3 1158 1626 (N32E11 L=183)  1094    205        2 (1102/1527-1211/1704) -39.8 E17-E04 N26-N36
      4  584 1680 (N20E51 L=143)   417    101        1 ( 554/1650- 614/1699)  -8.1 E54-E47 N17-N21
      5  808 1548 (N17E28 L=166)   554    166        1 ( 764/1504- 872/1609) -24.0 E32-E23 N15-N20
      6  739 1489 (N12E30 L=164)   449    145        0 ( 695/1441- 796/1541) -16.7 E34-E24 N10-N14
      7  424 1320 (S05E45 L=149)   767    228        0 ( 390/1217- 451/1416) -79.7 E48-E42 S11-N00
      8 1577 1285 (N24W26 L=220)   850    244        0 (1535/1198-1599/1391) -76.0 W21-W28 N17-N29
      9 1748 1273 (N30W41 L=235)   712    137        1 (1721/1194-1767/1326) -60.3 W37-W43 N25-N32
     10 1026 1266 (N06E07 L=187)   585    173        0 ( 976/1213-1105/1304)  54.5 E08-E02 N02-N09
     11  524 1230 (S08E35 L=159)   387    122        0 ( 484/1182- 549/1286)  81.9 E37-E33 S11-S04
     12 1587 1102 (N14W31 L=224)   225     58        0 (1568/1079-1606/1126)  90.0 W29-W32 N12-N15
     13 1063  909 (S12W05 L=199)   195     26        0 (1044/ 887-1083/ 934) -45.0 W04-W06 S13-S10
     14  189  862 (S37E63 L=131)   840    127        0 ( 142/ 828- 240/ 905)  17.5 E73-E54 S39-S33
     15 1003  666 (S27W10 L=203)   318     78        2 ( 967/ 646-1037/ 686)   9.5 W06-W12 S28-S25
     16 1478  343 (S29W59 L=253)   892     95        0 (1429/ 325-1534/ 369)  15.3 W54-W65 S30-S27
  

    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