Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1034 1747 (N51W29 L=290)   553     74        1 (1008/1712-1060/1780)   6.3 W25-W34 N48-N53
      2 1448 1717 (N35W58 L=319)   397     41        0 (1425/1702-1476/1735) -56.3 W55-W60 N33-N36
      3  918 1679 (N50W15 L=275)   365     48        0 ( 890/1658- 949/1697) -53.1 W12-W16 N48-N52
      4 1688 1612 (N20W73 L=333)  1358     80        3 (1654/1584-1728/1643) -33.7 W68-W79 N18-N22
      5  515 1357 (N39E28 L=233)  1452    290        4 ( 404/1280- 628/1392) -28.8 E39-E19 N31-N43
      6 1038 1328 (N24W09 L=269)  1074    349        1 ( 905/1300-1208/1356) -20.2 W00-W19 N20-N27
      7 1169 1233 (N15W14 L=274)   196     19        0 (1147/1218-1192/1252) -45.0 W12-W14 N14-N16
      8 1761 1150 (S06W51 L=312)  1408    161        5 (1723/1072-1789/1228)  45.0 W46-W56 S09-S01
      9  554 1076 (N20E28 L=233)   933    238        1 ( 467/ 998- 648/1137) -47.7 E34-E23 N14-N25
     10  205  958 (N19E62 L=199)  2623    439        5 (  98/ 811- 342/1096) -28.5 E85-E49 N08-N27
     11 1659  620 (S35W34 L=295)  6658    815       13 (1458/ 384-1849/ 959)  33.4 W14-W56 S49-S18
     12 1272  896 (S06W11 L=272)   286     64        0 (1251/ 860-1299/ 935) -90.0 W10-W12 S08-S03
     13 1216  857 (S07W07 L=268)   272    103        0 (1170/ 846-1274/ 868) -26.6 W04-W10 S08-S05
     14 1246  663 (S19W04 L=265)   646    268        1 (1180/ 564-1275/ 770)  48.0 E02-W07 S23-S12
     15 1340  312 (S49W02 L=262)   391     61        1 (1321/ 287-1362/ 347)  59.0 W00-W03 S51-S45
     16 1254  260 (S51E08 L=253)   674     79        1 (1219/ 233-1293/ 285) -38.7 E10-E04 S53-S47
     17  929  237 (S42E34 L=227)   498     78        0 ( 892/ 219- 970/ 260)   0.0 E37-E29 S42-S40
     18 1144  209 (S52E20 L=240)   536     76        1 (1111/ 186-1185/ 231)   0.0 E24-E15 S53-S50
  

    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