Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  587 1666 (N31E49 L=  8)   823    185        2 ( 572/1568- 612/1742) -45.0 E55-E42 N24-N36
      2  880 1487 (N25E19 L= 38)   230     54        0 ( 857/1463- 908/1506) -26.6 E20-E16 N23-N25
      3 1320 1435 (N30W12 L= 69)   438    212        0 (1227/1402-1403/1499)  -8.1 W04-W18 N28-N32
      4  762 1445 (N20E25 L= 32)   562    105        2 ( 717/1420- 817/1477)  33.7 E28-E21 N17-N21
      5  193 1225 (S04E67 L=350)   715     92        3 ( 164/1175- 221/1263) -24.4 E73-E62 S06-S01
      6 1698 1158 (N20W44 L=101)   607     81        3 (1655/1135-1739/1186)   0.0 W39-W47 N17-N21
      7  416  858 (S22E38 L= 19)  3157    661        5 ( 200/ 702- 642/1086)  36.7 E67-E24 S36-S04
      8 1751 1064 (N15W49 L=106)   440     76        1 (1713/1041-1779/1083)  23.2 W45-W52 N13-N16
      9  186 1008 (S17E63 L=354)  1067    116        2 ( 152/ 951- 209/1058) -63.4 E68-E59 S19-S13
     10  143  944 (S21E71 L=346)  2288    252        6 ( 111/ 880- 190/1043) -34.7 E80-E62 S25-S15
     11 1701  937 (N06W47 L=104)   317     86        0 (1675/ 899-1732/ 974) -29.7 W43-W50 N04-N08
     12 1508  884 (S00W33 L= 90)   553    104        1 (1452/ 857-1553/ 909)  15.9 W29-W36 S02-N00
     13 1505  741 (S09W36 L= 93)   293    163        0 (1476/ 715-1557/ 795)   0.0 W33-W40 S10-S05
     14 1241  707 (S16W19 L= 77)   383    142        1 (1181/ 686-1313/ 723)  35.5 W16-W23 S18-S13
     15 1705  578 (S15W60 L=117)  1755    282        2 (1614/ 494-1797/ 687)  59.0 W54-W65 S21-S06
     16  676  534 (S39E15 L= 42)   292     56        0 ( 659/ 502- 691/ 560) -56.3 E16-E13 S40-S36
     17 1479  393 (S31W51 L=108)   917    107        1 (1431/ 367-1528/ 432)  33.7 W47-W54 S33-S28
     18 1286  373 (S37W36 L= 93)   284     49        1 (1261/ 354-1309/ 392) -11.3 W33-W38 S37-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