Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  855 1720 (N38E33 L=142)   447     85        0 ( 822/1692- 889/1760) -24.0 E37-E28 N36-N40
      2  914 1576 (N29E21 L=154)   713    182        1 ( 817/1535- 976/1620)  52.1 E25-E17 N24-N33
      3 1659 1275 (N27W38 L=213)  2716    765        1 (1446/ 984-1793/1554) -39.6 W18-W53 N12-N41
      4  590 1437 (N13E37 L=137)   696    143        1 ( 544/1387- 641/1496)  76.0 E39-E35 N09-N17
      5 1320 1358 (N24W11 L=186)   496    228        0 (1235/1324-1396/1427)   0.0 W04-W16 N22-N26
      6 1159 1371 (N21W01 L=175)   245    117        0 (1115/1330-1200/1413) -33.7 E03-W03 N19-N23
      7  584 1339 (N07E34 L=140)   243     48        0 ( 557/1326- 617/1358)  11.3 E36-E31 N06-N08
      8 1362 1296 (N21W15 L=189)   230     62        0 (1340/1274-1383/1323)  71.6 W13-W16 N20-N23
      9  262 1229 (S05E58 L=117)   447     41        1 ( 243/1205- 284/1258)  90.0 E59-E55 S06-S03
     10  497 1025 (S13E33 L=141)   852    185        3 ( 445/ 931- 560/1071)  66.0 E37-E28 S18-S09
     11  254  846 (S29E54 L=120)   702    205        0 ( 214/ 785- 288/ 912)  48.4 E60-E49 S33-S24
     12  994  873 (S11W01 L=176)   394    135        0 ( 928/ 858-1050/ 888)   9.5 E02-W04 S12-S09
     13  716  806 (S21E14 L=160)   290     94        0 ( 677/ 781- 751/ 842)  51.3 E16-E12 S23-S18
     14 1676  732 (S05W50 L=224)   302     70        1 (1649/ 707-1717/ 748)  26.6 W47-W53 S06-S02
     15 1529  672 (S11W41 L=215)   364     63        1 (1499/ 642-1559/ 698) -26.6 W38-W42 S12-S09
     16 1628  579 (S14W53 L=227)  2226    336        7 (1579/ 465-1694/ 683) -60.3 W46-W59 S22-S08
     17 1097  607 (S25W14 L=189)   256     58        0 (1070/ 589-1131/ 627)  56.3 W12-W16 S26-S23
     18 1501  560 (S19W43 L=218)   615    168        1 (1425/ 535-1578/ 589)  31.0 W38-W48 S21-S15
     19  463  549 (S45E33 L=141)  1093    200        2 ( 396/ 504- 553/ 586)  27.8 E43-E24 S50-S40
  

    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