Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1275 1625 (N35E00 L=305)   867    159        2 (1213/1577-1349/1663) -11.3 E05-W05 N32-N36
      2  916 1446 (N14E16 L=289)   231     52        0 ( 885/1431- 949/1461)  14.0 E18-E14 N13-N14
      3  404 1412 (N02E49 L=257)   317     57        1 ( 388/1382- 416/1446) -63.4 E50-E46 N00-N04
      4  307 1363 (S02E56 L=250)   350     68        1 ( 278/1336- 329/1395)  53.1 E57-E54 S04-N00
      5 1560 1310 (N23W25 L=331)   270     68        0 (1548/1273-1574/1351) -68.2 W23-W26 N20-N25
      6 1811 1327 (N33W48 L=354)   394     41        0 (1781/1312-1839/1341)  14.0 W44-W52 N31-N34
      7 1428 1303 (N18W17 L=322)   410    103        0 (1379/1275-1482/1323)   8.1 W13-W20 N17-N19
      8 1812 1098 (N18W50 L=355)  1543    284        2 (1776/ 983-1845/1246) -68.2 W45-W54 N11-N27
      9  131 1088 (S20E70 L=236)  1962    138        3 (  97/1019- 160/1140)  66.8 E76-E65 S24-S17
     10 1770  946 (N08W49 L=354)   492     78        0 (1745/ 911-1804/ 982) -18.4 W46-W52 N07-N10
     11 1604  951 (N04W36 L=341)   273     85        0 (1561/ 937-1647/ 970)  38.7 W33-W38 N01-N05
     12  690  812 (S27E15 L=290)   501    122        0 ( 628/ 775- 741/ 846)  -7.1 E19-E11 S27-S25
     13  351  626 (S46E44 L=261)  1406    218        1 ( 254/ 586- 448/ 667)  -2.4 E57-E32 S48-S43
     14  557  575 (S45E22 L=284)   375     65        0 ( 521/ 546- 590/ 599) -15.9 E25-E18 S46-S43
     15 1795  512 (S12W75 L= 21)  1710    103        3 (1764/ 463-1836/ 552) -71.6 W69-W82 S14-S08
     16  842  407 (S49W07 L=312)   253     50        0 ( 812/ 385- 870/ 432)  -9.5 W03-W09 S49-S47
     17 1583  398 (S25W61 L=  6)   659     83        1 (1536/ 372-1616/ 420)   5.7 W55-W65 S26-S23
     18 1449  247 (S38W67 L= 12)  1223    121        0 (1384/ 223-1511/ 272)  13.2 W57-W76 S39-S35
  

    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