Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  443 1621 (N13E58 L= 18)  1576    388        6 ( 392/1579- 515/1685) -15.5 E68-E50 N10-N17
      2  858 1590 (N21E27 L= 48)   308    113        3 ( 808/1552- 892/1637) -26.6 E31-E23 N19-N23
      3  374 1500 (N05E57 L= 18)  1008    202        4 ( 339/1442- 408/1571) -71.6 E60-E52 N00-N09
      4 1292 1354 (N19W06 L= 81)   908    262        1 (1267/1230-1319/1478) -68.2 W02-W10 N12-N27
      5 1626 1345 (N30W28 L=104)   307     32        0 (1605/1328-1651/1366)   0.0 W26-W30 N28-N31
      6  251 1320 (S08E61 L= 14)   447     90       -1 ( 235/1290- 269/1349) -71.6 E63-E58 S09-S06
      7 1482 1188 (N16W22 L= 97)   367    100        1 (1464/1134-1500/1241) -56.3 W19-W23 N12-N18
      8  188  952 (S31E57 L= 18)  1715    412        3 ( 100/ 923- 286/ 982)  14.7 E66-E47 S33-S27
      9 1303  912 (S05W18 L= 94)   229     46        2 (1280/ 881-1324/ 939) -76.0 W16-W19 S06-S02
     10  734  867 (S23E13 L= 63)   198     60        2 ( 706/ 845- 757/ 895)  76.0 E13-E11 S24-S20
     11 1855  862 (N11W60 L=136)   449     73        0 (1843/ 824-1869/ 894) -36.9 W57-W62 N08-N12
     12  465  741 (S38E29 L= 46)  1466    334        3 ( 347/ 709- 586/ 792)  18.4 E40-E19 S42-S33
     13 1834  666 (N00W69 L=144)  1232    241        2 (1797/ 642-1887/ 704)  14.9 W63-W79 S01-N03
     14  965  593 (S32W10 L= 85)   986    199        0 ( 882/ 522-1050/ 680)  74.1 W06-W12 S38-S24
     15 1302  579 (S22W30 L=106)   551    155        1 (1222/ 543-1367/ 638)  66.0 W26-W33 S26-S17
     16 1569  579 (S14W48 L=123)   557    180        2 (1504/ 549-1626/ 619)   5.2 W41-W52 S15-S11
     17 1282  431 (S31W37 L=112)   673    164        1 (1194/ 411-1359/ 465)  52.1 W31-W40 S35-S26
  

    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