Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1124 1796 (N52W42 L=289)   595     87        1 (1087/1760-1147/1833)  24.4 W37-W48 N48-N53
      2  663 1421 (N40E14 L=233)   982    199        2 ( 586/1368- 755/1462) -26.6 E20-E08 N33-N42
      3 1212 1415 (N25W22 L=269)  1037    303        3 (1085/1387-1367/1447) -18.4 W14-W32 N20-N27
      4  542 1412 (N42E25 L=222)   262     46        1 ( 524/1386- 563/1436) -76.0 E26-E23 N40-N44
      5 1083 1303 (N22W11 L=258)   303    120        1 (1050/1266-1131/1349) -68.2 W09-W12 N18-N25
      6  978 1223 (N20W02 L=249)   230    100        0 ( 929/1202-1025/1252)   0.0 E01-W05 N18-N20
      7 1854 1189 (S07W63 L=310)   452     37        1 (1833/1161-1871/1211)  33.7 W60-W66 S08-S05
      8  729 1144 (N21E15 L=232)   817    240        1 ( 631/1069- 821/1200) -47.7 E20-E10 N14-N25
      9  288 1015 (N22E52 L=195)   827    100        2 ( 243/ 975- 329/1052) -36.9 E56-E47 N18-N24
     10 1888  914 (S24W62 L=309)  2086    212        3 (1841/ 817-1928/1011)  39.8 W54-W70 S29-S19
     11 1402  675 (S23W14 L=261)   362    167        0 (1338/ 628-1453/ 735)  24.0 W09-W18 S24-S20
     12 1724  566 (S43W43 L=290)  3257    396        5 (1594/ 443-1857/ 706)  20.3 W27-W62 S49-S35
     13  516  502 (S13E47 L=200)   509    124        0 ( 478/ 449- 550/ 552) -90.0 E48-E45 S16-S09
     14 1452  350 (S50W13 L=260)   395     53        0 (1436/ 331-1473/ 377) -80.5 W11-W15 S52-S46
     15 1382  320 (S50W05 L=252)   500     53        0 (1360/ 296-1414/ 340)   9.5 W02-W08 S51-S47
     16 1079  283 (S42E20 L=227)   654     89        1 (1039/ 256-1119/ 323)  14.0 E23-E15 S44-S39
     17  599  276 (S30E58 L=189)   680     72        0 ( 553/ 263- 631/ 293)  -6.3 E63-E53 S31-S29
     18 1257  249 (S52E08 L=239)   525     68        1 (1226/ 222-1294/ 272)  14.0 E12-E04 S54-S48
  

    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