Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1405 1734 (N40W54 L=260)   542     55        1 (1369/1713-1436/1751)  -8.1 W50-W58 N38-N42
      2 1222 1712 (N45W37 L=243)   276     49        0 (1200/1691-1249/1728)   9.5 W33-W39 N43-N46
      3 1017 1677 (N48W18 L=224)   303     67        0 ( 995/1642-1045/1702)  33.7 W15-W21 N45-N49
      4  699 1554 (N47E13 L=193)  1487    319        6 ( 593/1490- 856/1616)   5.0 E22-W01 N43-N49
      5 1574 1369 (N15W43 L=250)   582    212        0 (1488/1298-1627/1435) -65.6 W39-W45 N09-N20
      6  510 1347 (N36E30 L=176)   321     78        0 ( 469/1324- 545/1368) -29.7 E34-E27 N33-N37
      7  897 1145 (N16E05 L=201)  1069    355        2 ( 821/1009- 985/1278) -74.5 E07-E02 N07-N25
      8  606 1232 (N27E23 L=183)   302     69        0 ( 570/1212- 654/1245)  -9.5 E25-E19 N25-N27
      9 1583  882 (S13W32 L=238)   230     67        0 (1557/ 852-1600/ 921)  53.1 W29-W33 S14-S10
     10  311  773 (N03E53 L=153)   420     79        1 ( 284/ 735- 338/ 803) -68.2 E55-E50 N00-N05
     11 1651  727 (S25W36 L=242)  1024    200        1 (1581/ 687-1729/ 784) -24.8 W30-W43 S27-S19
     12  721  736 (S04E24 L=183)   245     71        0 ( 710/ 700- 732/ 778)  68.2 E24-E22 S06-S01
     13 1773  535 (S44W55 L=261)   583     44        2 (1744/ 519-1800/ 555) -19.7 W48-W62 S46-S40
     14 1604  441 (S45W34 L=240)  2390    390        4 (1435/ 347-1748/ 533)   4.9 W16-W51 S49-S41
     15  351  442 (S19E67 L=140)  3314    224        8 ( 289/ 375- 431/ 518) -63.4 E75-E60 S24-S13
     16  608  417 (S23E42 L=164)  1026    194        3 ( 580/ 327- 639/ 501)  47.5 E48-E36 S29-S17
     17  470  476 (S17E50 L=156)   349     52        1 ( 447/ 448- 491/ 499) -90.0 E52-E48 S19-S15
     18 1047  368 (S35E14 L=193)   434    126        2 (1026/ 307-1074/ 423)  68.2 E16-E10 S40-S30
  

    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