Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1221 1664 (N41W33 L=254)  1588    367        2 (1051/1625-1365/1701) -10.5 W18-W45 N38-N44
      2  623 1535 (N47E20 L=201)  2396    598        5 ( 416/1448- 910/1648)   4.8 E42-W07 N43-N48
      3 1787 1473 (N14W70 L=292)   811     99        0 (1746/1440-1834/1505) -29.1 W66-W76 N11-N16
      4 1387 1285 (N15W28 L=249)   441    226        0 (1303/1201-1440/1355) -71.6 W24-W29 N08-N20
      5  368 1240 (N31E43 L=178)   290     83        0 ( 358/1198- 383/1283)  80.5 E44-E42 N27-N33
      6  635 1169 (N22E22 L=200)   339     84        1 ( 612/1130- 660/1210) -80.5 E23-E20 N19-N25
      7 1772 1076 (S06W51 L=272)   408     77        1 (1731/1053-1810/1105) -35.5 W47-W54 S08-S03
      8  703  995 (N11E20 L=202)   462    113        2 ( 664/ 942- 738/1043) -74.1 E21-E17 N07-N14
      9  361  975 (N14E45 L=177)   349     78        1 ( 337/ 937- 391/1000)  31.0 E47-E42 N12-N15
     10 1361  783 (S14W16 L=237)   280    102        0 (1343/ 738-1385/ 832)  51.3 W13-W17 S16-S11
     11  502  688 (S04E39 L=182)   410    114        1 ( 490/ 630- 521/ 747)  76.0 E41-E37 S08-N00
     12 1427  650 (S24W18 L=240)   605    142        0 (1389/ 614-1476/ 707) -54.5 W15-W21 S27-S19
     13  771  544 (S17E25 L=196)   298    106        0 ( 735/ 498- 791/ 586)  39.8 E29-E23 S19-S14
     14 1585  449 (S44W32 L=253)  3202    476        3 (1402/ 354-1780/ 574)   5.6 W13-W55 S47-S39
     15 1022  385 (S33E14 L=207)   286     88        0 ( 989/ 353-1050/ 426)  39.8 E17-E11 S35-S30
     16  452  384 (S24E59 L=163)   952     78        3 ( 430/ 338- 475/ 421)  24.8 E66-E53 S27-S21
     17 1317  308 (S48W05 L=227)   393     70        0 (1280/ 296-1353/ 321) -26.6 W02-W08 S49-S45
     18  959  149 (S54E36 L=185)   492     44        0 ( 934/ 134- 987/ 165) -45.0 E39-E33 S56-S51
  

    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