Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1706 1343 (N06W53 L=259)   303     71        0 (1693/1305-1724/1373)  63.4 W51-W56 N03-N07
      2  317 1229 (N35E47 L=158)   825    207        3 ( 233/1185- 410/1272) -18.4 E58-E37 N30-N37
      3 1093 1148 (N12W07 L=213)   183     41        2 (1069/1125-1112/1167) -71.6 W05-W08 N10-N13
      4 1495 1093 (S01W29 L=235)   502    115        2 (1429/1067-1558/1112) -12.5 W24-W33 S02-N00
      5 1633  831 (S20W33 L=239)   806    286        1 (1530/ 726-1723/ 903)  12.5 W23-W41 S24-S17
      6 1089  801 (S07E02 L=204)   369    149        2 (1048/ 741-1136/ 858) -90.0 E03-W00 S11-S03
      7 1676  740 (S28W36 L=241)   702    208        0 (1636/ 651-1706/ 830)  51.3 W31-W39 S32-S22
      8 1089  658 (S15E05 L=200)   249     66        1 (1064/ 625-1120/ 685)  76.0 E06-E03 S17-S13
      9 1345  598 (S26W09 L=214)   268     96        1 (1303/ 572-1381/ 637)   9.5 W05-W11 S27-S24
     10  978  504 (S22E16 L=189)   918    288        3 ( 869/ 434-1062/ 601)  20.6 E25-E09 S25-S17
     11 1168  474 (S29E06 L=200)   231    108        0 (1147/ 434-1203/ 514)  38.7 E08-E02 S31-S27
     12  486  411 (S18E55 L=150)   978    120        2 ( 459/ 345- 513/ 453)  48.8 E60-E51 S22-S14
     13 1363  405 (S41W05 L=211)   359     92        1 (1328/ 376-1403/ 435)   8.1 W02-W09 S42-S39
     14  447  341 (S22E67 L=139)   848     65        2 ( 418/ 323- 481/ 368) -18.4 E71-E62 S24-S20
     15  728  326 (S29E41 L=164)   495     63        2 ( 706/ 299- 755/ 354)  18.4 E44-E38 S30-S26
     16 1128  249 (S46E19 L=187)   390     79        0 (1096/ 234-1166/ 271)   8.1 E21-E14 S47-S44
  

    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