Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1220 1781 (N49W21 L=282)  1356    178        5 (1154/1731-1303/1840) -36.0 W16-W27 N43-N54
      2 1615 1661 (N38W55 L=315)  1297    194        1 (1516/1638-1698/1682)   0.0 W44-W72 N36-N39
      3 1541 1482 (N25W38 L=298)  7963   1115       13 (1113/1178-1901/1649) -20.7 W07-W69 N05-N37
      4 1531 1363 (N17W35 L=295)   420    106        1 (1515/1308-1553/1414) -90.0 W33-W36 N13-N20
      5  824 1348 (N18E12 L=248)   557    120        3 ( 774/1306- 876/1385) -45.0 E14-E08 N14-N20
      6 1307 1291 (N13W18 L=278)   291     77        0 (1271/1252-1347/1322) -45.0 W16-W20 N10-N14
      7 1747 1250 (N10W51 L=311)   429     75        2 (1734/1203-1774/1291)  80.5 W49-W53 N06-N12
      8 1660 1257 (N10W43 L=303)   324     50        0 (1636/1234-1681/1287)  36.9 W41-W45 N08-N12
      9  411 1206 (N11E40 L=220)   385    145        0 ( 359/1179- 474/1242)  20.6 E43-E35 N09-N12
     10  978 1193 (N08E02 L=258)   271     81        0 ( 943/1160-1012/1227) -45.0 E04-W00 N05-N09
     11 1894  873 (S13W67 L=327)  3507    279        2 (1845/ 764-1952/1018)  55.5 W62-W75 S19-S03
     12 1606  788 (S18W39 L=299)   995    321        0 (1551/ 648-1678/ 922)  76.8 W35-W43 S26-S09
     13  411  768 (S16E42 L=218)   822    215        3 ( 313/ 736- 502/ 791)   0.0 E50-E34 S17-S13
     14  397  602 (S26E48 L=212)   531     90        1 ( 366/ 571- 439/ 639)  24.0 E52-E43 S27-S23
     15 1272  589 (S30W16 L=276)   242     84        0 (1245/ 553-1308/ 621)  45.0 W14-W18 S32-S28
     16 1736  506 (S36W63 L=323)  1671    205        0 (1648/ 472-1808/ 559)  15.5 W50-W73 S38-S32
     17 1472  465 (S39W34 L=295)   653     97        2 (1424/ 425-1511/ 505) -18.4 W29-W38 S42-S36
     18 1156  449 (S40W08 L=268)   832    172        0 (1073/ 420-1244/ 486)  13.0 W01-W14 S41-S37
     19  373  383 (S40E69 L=191)   832     75        1 ( 338/ 353- 402/ 418)   7.1 E83-E59 S42-S38
  

    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